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Abstract 

Catalytic  self-ignition  of  methanol  over  platinum  supported  on  anodized  alumina  in  a  microburner  is  demonstrated.  The  effects  of  air  flow  rate 
and  air/methanol  ratio  on  the  system  start-up  time  and  steady  state  temperature  have  been  studied.  The  air  flow  rate  strongly  affects  the  system 
start-up  time.  However,  at  the  same  flow  rate,  the  start-up  time  is  not  affected  by  the  air/methanol  ratio  over  the  range  studied.  Computational  fluid 
dynamics  simulations  confirm  these  experimental  findings.  It  is  proposed  that  high  flow  rates  and  relatively  lower  fuel  content  can  lead  to  attaining 
steady  state  faster  with  minimal  fuel  utilization.  Transient  axial  temperature  profiles  showed  that  the  hot  spot  of  the  reaction  started  at  the  inlet  of 
the  burner  and  did  not  shift  downstream  with  time  under  most  reaction  conditions.  The  burner  effective  thermal  conductivity  was  increased  using 
copper  thermal  spreaders  on  the  burner  outer  walls,  allowing  for  thermally  uniform  walls.  The  uniform  temperature  profile  allowed  for  optimal 
integration  of  the  microburner  with  a  thermoelectric  device.  The  integrated  microburner/thermoelectric  device  was  shown  to  self- start  from  room 
temperature.  The  maximum  power  generated  with  the  thermoelectric  was  0.65  W  and  the  maximum  thermal  efficiency  was  1.1%. 

©  2008  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

There  is  an  increasing  demand  for  portable  power  supply 
using  higher  energy  density  systems  and  for  longer  operational 
time  between  recharges  [1].  Unfortunately,  battery  technol¬ 
ogy  has  not  kept  pace  with  the  increasing  power  demand  for 
portable  electronics,  which  forced  manufacturers  to  introduce 
power-down  options  to  prolong  the  battery  operational  time. 
Similarly,  the  military  has  also  been  interested  in  portable  power 
sources  for  long-duration  missions  where  portable  electronics 
are  increasingly  being  used  [2].  Hydrocarbons  possess  high 
gravimetric  and  volumetric  energy  density  compared  to  batter¬ 
ies.  For  instance,  the  energy  density  of  diesel  fuel  (12  kW  h  kg- 1 ) 
and  methanol  (5.6  kW  h  kg-1 )  are  each  at  least  an  order  of  mag¬ 
nitude  greater  than  that  of  lithium-ion  batteries  (0.2  kW  h  kg-1) 
[3,4].  In  recent  years,  new  approaches,  ranging  from  direct 
methanol  fuel  cells  [5,6],  to  reformed  hydrogen  fuel  cells 
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[2-4,7],  to  (thermal  to  electric  energy)  devices  integrated  with 
microburners  [8-12]  have  been  investigated  as  alternative  tech¬ 
nologies  that  generate  electric  power  from  fuels. 

Proton  exchange  membrane  (PEM)  fuel  cell  systems  operat¬ 
ing  on  pure  hydrogen  possess  high  energy  density  and  exhibit 
efficiencies  as  high  as  70%.  Hydrogen  production  from  fuels 
currently  appears  as  the  only  viable  option  due  to  the  low 
volumetric  energy  density  of  hydrogen  and  the  difficulty  and 
safety  issues  associated  with  its  storage.  Unfortunately,  the  inte¬ 
gration  of  a  fuel  reformer,  combustor  and  carbon  monoxide 
conversion  and  cleaning  units  (shift  reactors  plus  preferential 
oxidation  reactor  or  a  separation  unit)  with  a  fuel  cell  severely 
degrades  the  overall  system  performance  and  make  a  portable 
system  quite  complicated  (due  to  multiple  fuel  streams  and 
different  operating  temperatures  of  various  systems).  Holladay 
et  al.  reported  obtaining  lOOmW  maximum  electrical  output 
from  a  mesoscale  fuel  cell  operating  on  pure  hydrogen.  Based 
on  the  reported  H2  flow  rate  and  the  lower  heating  value  of 
H2,  we  calculated  the  thermal  efficiency  of  the  fuel  cell  to  be 
about  33%.  However,  the  efficiency  of  the  fuel  cell  dropped  to 
16.3%  when  integrated  with  the  methanol  reformer  and  com- 
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Fig.  1.  Methanol  delivery  system. 

bustor  and  the  maximum  electrical  power  produced  by  the 
integrated  system  was  23  mW.  The  authors  attributed  the  loss 
in  performance  to  the  operation  with  reformate  (H2,  CO2  and 
CO)  instead  of  pure  H2  [4].  The  highest  thermal  efficiency 
reported  for  their  integrated  reformer/combustor  without  the  fuel 
cell  was  9%  [3].  Therefore,  we  estimate  the  total  efficiency  of 
the  integrated  reformer/combustor/fuel  cell  system  to  be  about 
16.3%  x  9%  =  1.5%.  A  thermoelectric  device  integrated  with  a 
microburner  clearly  offers  a  simpler  design  over  a  reformed  fuel 
cell  system  at  a  comparable  efficiency,  as  was  shown  by  Federici 
et  al.  [9]. 

In  our  earlier  work,  we  found  that  hydrogen  self-ignites 
on  platinum  (Pt)  nanoparticles  dispersed  on  anodized  alumina 
wafers.  This  feature  simplifies  start-up  since  no  other  ignition 
sources  need  to  be  integrated.  Hydrogen  can  also  be  used  to 
ignite  hydrocarbon  fuels  [13].  However,  self-ignition  of  small 
hydrocarbons  (alkanes  up  to  butane)  has  not  been  feasible  due 
to  the  difficulty  in  activating  the  C-H  bond  cleavage.  Methanol 
has  been  shown  to  self-ignite  catalytically  on  supported  Pt  and 
Pd  catalysts  [14-16].  Self-ignition  of  methanol  has  also  been 
shown  to  occur  over  Pt/ALCh  even  when  the  burner  was  kept  at 
0°C  [16]. 

In  this  work,  we  study  the  catalytic  ignition  of  lean 
methanol/air  mixtures  over  Pt  catalyst  dispersed  on  anodized 
alumina  wafers  embedded  in  a  microbumer.  Our  focus  on  the 
simplest  alcohol  has  in  part  been  motivated  from  the  recent 
interest  in  oxygenated  fuels,  some  of  which  could  be  produced 
from  bio-renewable  sources,  and  its  inherent  advantage  in  terms 
of  being  liquid  at  room  temperature  but  with  a  relatively  high 
vapor  pressure.  The  effects  of  air  flow  rate  and  air/methanol 
ratio  on  the  start-up  time  and  steady  state  temperature  have  been 
studied  in  order  to  develop  strategies  for  optimal  start-up.  Com¬ 
putational  fluid  dynamics  (CFD)  simulations  are  carried  out  to 
confirm  the  experimental  findings.  Finally,  control  of  temper¬ 
ature  profiles  is  demonstrated  and  the  efficiency  of  integrated 
self- starting,  microburner- thermoelectric  devices  is  estimated 
from  experimental  data. 

2.  Experimental 

2.7.  Microburner  setup 

Fig.  1  shows  the  methanol/air  mixing  system.  The  air  flows 
through  two  methanol  bubblers  where  it  gets  enriched  with 
methanol  vapor.  The  first  bubbler  was  held  at  room  temperature 
while  the  second  was  held  either  at  0  °C  or  at  room  temperature. 

Fig.  2  shows  a  schematic  of  the  microburner  design.  The 
burner  consists  of  two  rectangular  sheets  of  0.79  mm  thick  316- 
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Fig.  2.  Schematic  of  the  burner  design  (side  view). 

alloy  stainless  steel,  with  a  bolt  pattern  around  the  exterior.  A 
thin  (540  p,m)  gasket  separates  the  two  sheets  to  form  a  rect¬ 
angular  microchannel  approximately  6  cm  long  by  1  cm  wide. 
Stainless  steel  inlet  and  outlet  tubes  were  welded  to  one  of  the 
exterior  sheets.  Fine  metal  screens  were  spot  welded  near  the 
inlet  and  outlet  of  the  channel  in  order  to  induce  uniform  flow 
over  the  microchannel.  As  in  previous  work  [8],  two  copper 
thermal  spreaders  were  used  to  change  the  effective  thermal 
conductivity  of  the  walls  in  order  to  facilitate  axial  heat  transfer. 
In  order  to  provide  a  uniform  surface  of  the  same  dimensions 
of  the  thermoelectric,  a  secondary  copper,  thermal  spreader  was 
installed  in  experiments  conducted  with  thermoelectrics. 

An  AI2O3  wafer,  approximately  100  jim  in  thickness,  has 
been  inserted  into  the  burner,  which  reduced  the  burner  gap  size 
to  ~440  |mm.  Eight  K-type  thermocouples  were  placed  on  the 
outer  wall  of  the  burner  to  record  the  temperature  profile.  The 
temperature  profile  during  ignition  was  monitored  via  Lab  view® 
software  and  National  Instruments  data  acquisition  board.  The 
burner  setup  has  been  described  in  detail  in  previous  work 
[8,9]. 

2.2.  Catalyst  preparation 

The  AI2O3  support  was  prepared  by  completely  anodizing 
100  [Jim-thick  1145  alloy  aluminum  foil  to  create  thin  alumina 
wafers  that  serve  as  high-surface  area  porous  catalyst  supports 
(approximately  100  p,m  thick).  The  foil  is  immersed  in  0.3  M 
aqueous  oxalic  acid  at  a  potential  of  40  V  for  48  h.  The  system 
is  held  at  approximately  1-2  °C  using  an  ice  bath.  The  pore 
structure  is  semi-ordered  with  an  average  pore  diameter  of  50  nm 
and  a  surface  area  of  approximately  14  m2  g-1 .  The  pores  travel 
the  thickness  of  the  wafer. 

The  Pt/Al203  catalyst  was  prepared  by  wet  impregnation, 
where  the  platinum  was  deposited  by  immersing  the  AI2O3 
wafer  in  0.007  M  aqueous  dihydrogen  hexachloroplatinate(IV) 
with  pH  3  (for  optimal  surface  charge)  for  3  h  [17].  The  AI2O3 
wafers  were  then  removed  from  the  solution  and  dried  overnight, 
and  subsequently  calcined  in  air  at  500  °C  for  3  h.  Prior  to  reac¬ 
tivity  measurements,  the  catalyst  was  reduced  in  100  seem  10% 
H2  (balance  He)  at  300  °C  for  3h.  The  Pt  weight  loading  was 
approximately  1  wt.%  as  determined  by  atomic  absorption  spec¬ 
troscopy.  The  catalyst  loading  in  the  burner  was  between  75  and 
85  mg. 
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2.3.  Composition  measurements 

The  exhaust  gases  were  sampled  with  a  gas-sampling  valve, 
and  the  composition  was  monitored  with  a  HP  Series  6890  gas 
chromatograph  (GC),  using  both  a  thermal  conductivity  detector 
and  a  flame  ionization  detector.  As  seen  in  Fig.  1,  the  air  flows 
through  two  methanol  bubblers  where  the  incoming  stream  gets 
saturated  with  methanol  at  the  low  air  flow  rates  as  confirmed 
by  the  GC.  However,  at  the  higher  flow  rates  the  temperature  of 
the  bubbler  drops  due  to  the  fast  evaporation  of  methanol  and 
it  is  not  clear  whether  the  incoming  stream  is  saturated  or  not. 
The  methanol  mole  fraction  in  the  mixture  varied  from  3.6  to 
12.3%.  The  fraction  of  methanol  in  the  mixture  was  determined 
by  back  calculating  the  methanol  percentage  from  the  products 
composition  assuming  the  total  oxidation  reaction  stoichiometry 
(this  is  indeed  the  product  distribution  as  the  results  show  below). 
To  verify  this  method,  we  later  measured  the  methanol  fraction 
in  the  inlet  mixture  directly  and  the  composition  was  within  10% 
of  the  calculated. 

A  methanol  mole  fraction  of  12.3%  in  the  mixture  is  equiva¬ 
lent  to  stoichiometric  oxygen  to  methanol  ratio  of  1 .5.  Therefore, 
all  experiments  were  conducted  under  fuel  lean  conditions.  The 
methanol  conversion  was  calculated  based  on  a  carbon  bal¬ 
ance  using  the  measured  composition  of  carbon  dioxide  and 
methanol,  if  any,  in  the  exhaust  gases.  The  error  in  the  total 
mass  balance  was  lower  than  10%  in  all  experiments. 

3.  Results 

3. 1 .  Catalyst  characterization 

Scanning  transmission  electron  microscopy  (STEM)  was  per¬ 
formed  on  a  JEOL  201  OF  FASTEM  field  emission  gun  scanning 
transmission  electron  microscope  (STEM).  The  catalyst  wafer 
used  in  the  burner  integrated  with  the  thermoelectric  was  taken 
out  of  the  burner  for  characterization.  The  catalyst  wafer  was 
used  for  about  60  h  on  stream  with  over  20  ignition  and  shut 
down  cycles.  No  apparent  catalyst  deactivation  was  observed. 
No  deactivation  was  seen  in  previous  microreactor  work  involv¬ 
ing  longer  runs  and  more  adverse  temperature  conditions  than 
the  ones  used  herein  [8].  The  catalyst  wafer  was  crushed  into 
fine  powder  and  dispersed  in  ethanol.  A  sample  of  the  dispersed 
catalyst  was  dropped  on  a  carbon  coated  copper  TEM  grid. 

We  used  high  angle  annular  dark  field  imaging  (HAADF) 
in  the  STEM  mode  to  enable  detection  of  the  smallest  parti¬ 
cles.  Fig.  3  shows  a  representative  STEM  image  of  the  spent 
catalyst.  The  Pt  nanoparticles  and  the  AI2O3  pores  can  be  seen 
in  the  image  as  bright  and  dark  areas,  respectively.  Over  sev¬ 
eral  images,  we  measured  the  particle  size  of  more  than  900 
Pt  nanoparticles.  The  corresponding  particle  size  distribution 
is  shown  in  Fig.  4.  The  Pt  number  average  particle  diameter 
(Aavg  =  JA  n i  di  /  JA  n i ,  di  being  the  particle  diameter  and  ni 
the  number  of  particles)  was  4.5  nm.  However,  the  particle  size 
distribution  is  broad  with  particles  as  large  as  25  nm.  The  Pt 
dispersion  was  about  10%  based  on  the  mean  surface  diam¬ 
eter  (As  =  J2inidf  /Yhfiidf)  [18].  The  low  dispersion  is  due 
to  sintering  from  operating  at  temperatures  as  high  as  650  °C. 


Fig.  3.  STEM  image  of  the  spent  catalyst  (1  wt.%  Pt/AbCb). 


However,  despite  the  broad  particle  size  distribution  and  low 
dispersion,  ignition  at  room  temperature  was  achieved  under  all 
conditions  studied  even  after  the  catalyst  has  been  on  stream  for 
60  h. 

3.2.  Effect  of  flow  rate  and  methanol  mole  fraction  on 
temperature  profiles 

Figs.  5  and  6  show  the  transient  temperature  profiles  during 
methanol  ignition  at  1000  seem  air  flow  rate  and  10%  methanol 
mole  fraction  in  the  mixture.  The  hot  spot  reached  100  °C  within 
10  s  of  start-up.  The  transient  temperature  profile  at  position 
1  cm  (see  Fig.  5)  was  fitted  to  a  first-order  dynamics  model, 
a{  1  —  e_r/T),  where  a  is  a  normalizing  factor,  t  is  time  and  t  is  the 
time  constant  of  the  system.  Accordingly,  we  find  that  the  sys- 
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1000  seem,  10%  MeOH 


Fig.  5.  Transient  temperature  profiles.  Air  flow  rate  1000  seem,  10%  methanol 
in  the  mixture.  Thermocouple  positions  stated  with  the  catalyst  starting  at  0  cm. 


tem  reached  95%  of  the  steady  state  temperature  (at  time  =  3r) 
in  about  9  min.  The  location  of  the  hot  spot  near  the  entrance 
of  the  catalytic  region  indicates  very  fast  kinetics  of  methanol 
oxidation  on  Pt.  The  reaction  products  were  analyzed  at  both 
approximately  60  s  from  starting  the  experiment  and  at  steady 
state.  Both  conditions  showed  100%  methanol  conversion  to 
water  and  carbon  dioxide. 

Fig.  7  shows  the  steady  state  temperature  profiles  at  differ¬ 
ent  methanol  mole  fractions  and  air  flow  rates.  At  low  methanol 
mole  fractions,  the  temperature  profiles  are  more  uniform  and 
the  average  temperature  is  much  lower  than  at  high  methanol 
mole  fractions.  The  hot  spot  was  near  the  entrance  (about  5  mm) 
of  the  catalytic  region  for  all  conditions  with  the  exception  of 
1000  seem  air  flow  rate  where  the  hot  spot  shifted  slightly  down¬ 
stream.  At  all  flow  rates,  the  methanol  conversion  was  >99% 
and  the  only  products  detected  at  steady  state  were  water  and 
carbon  dioxide.  During  the  first  minute  of  starting  the  experi¬ 
ment,  formaldehyde  and  methanol  breakthrough  were  observed 
at  low  methanol  mole  fractions.  No  methanol  breakthrough  was 


Thermocouple  position  (cm) 

Fig.  6.  Transient  temperature  profiles.  Air  flow  rate  =  1000  seem,  10%  methanol 
the  mixture. 


Thermocouple  Position  (cm) 

Fig.  7.  Effect  of  air  flow  rate  and  methanol  mole  fraction  in  the  mixture  on 
steady  state  temperature  profiles. 

detected  at  60  s  at  higher  methanol  mole  fractions.  This  is  most 
probably  due  to  the  higher  temperatures  reached  with  higher 
methanol  mole  fraction  in  the  feed. 

3.3.  Effect  of  flow  rate  and  methanol  mole  fraction  on 
start-up 

The  steady  state  temperatures  vary  with  varying  flow  rate  and 
methanol  composition.  In  order  to  compare  the  system  start¬ 
up  time  for  different  conditions,  the  dimensionless  temperature 
(T/Tss)  was  introduced,  as  the  temperature  T  at  position  of  1  cm 
divided  by  the  steady  state  temperature,  Tss,  at  the  same  position. 

At  the  same  air  flow  rate,  higher  methanol  mole  fraction 
resulted  in  a  faster  heating  rate;  however,  the  steady  state  tem¬ 
perature  was  also  higher.  Therefore,  as  shown  in  Fig.  8,  the  time 
it  takes  the  system  to  reach  steady  state  (start-up  time)  is  weakly 
affected  by  the  air/methanol  ratio.  The  air  flow  rate  on  the  other 
hand,  highly  affects  the  start-up  time,  as  shown  in  Fig.  8.  At 
higher  air  flow  rate,  the  rate  of  heat  release  by  methanol  oxi¬ 
dation  increases,  and  therefore,  the  temperature  increases  at  a 
faster  rate.  Jiang  et  al.  [16]  have  reported  similar  results  on  the 
effect  of  flow  rate  at  a  methanol  mole  fraction  of  1 1.5%. 

The  above  results  can  be  rationalized  by  simplifying  the  prob¬ 
lem  to  a  continuous  stirred  tank  reactor  and  solving  the  time 


Time  (sec) 

Fig.  8.  Effect  of  methanol  mole  fraction  in  the  mixture  on  system  start-up  in 
terms  of  dimensionless  temperature  T/Tss  at  1  cm  (from  the  catalyst  entrance) 
vs.  time. 
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dependent  energy  balance  for  an  ideal  gas  (assuming  complete 
fuel  conversion).  Eq.  (1)  gives  the  temperature  versus  time  in 
terms  of  the  total  flow  rate  Vm,  initial  temperature  7),  methanol 
mole  fraction  yMeOH,  mass  of  the  system  m  (this  is  mainly 
comprised  of  the  solid  walls  in  a  microchemical  system),  heat 
capacity  of  the  system  Cp,  heat  capacity  of  the  fluid,  Cpf,  den¬ 
sity  of  the  fluid  at  the  inlet,  pm,  heat  of  reaction  A77,  heat  loss 
coefficient  h,  inlet  methanol  concentration  Cm,  and  surface  area 
of  the  system  available  for  heat  transfer  A  (the  term  pinCpf  is 
assumed  to  be  independent  of  temperature,  a  common  assump¬ 
tion  in  the  combustion  literature;  a  reference  temperature  for  the 
gas  enthalpy  and  heat  of  reaction  of  298  K  is  chosen). 

^intPin^Pf^in  AH  yMeOH^in]  T"  Tdm\)hA 
hA  +  VjnPinCpf 

^intPinCpf^in  ~\~  AH  yMeOH^in]  T^^hA 
hA  V'm P'mCpf 

The  characteristic  time  scale  of  the  system  is  given  by 
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Fig.  9.  Steady  state  temperature  profiles  with  copper  thermal  spreaders  installed 
on  the  outer  walls  of  the  burner. 


hA  +  V(pCP)0Ut 

The  time  constant  of  the  system,  r,  does  not  depend  on  the 
methanol  mole  fraction.  However,  t  is  inversely  proportional 
to  the  total  flow  rate.  The  above  expression  for  r  confirms  the 
observed  experimental  results,  where  the  system  start-up  was 
only  affected  by  the  air  flow  rate  and  insensitive  to  the  methanol 
mole  fraction.  Two-dimensional  (in  space),  time  dependent 
CFD  simulations,  summarized  in  Appendix  A,  indicate  that  this 
dependence  on  flow  rate  and  independence  of  mixture  compo¬ 
sition  is  probably  a  generic  behavior  of  catalytic  combustion. 

3.4.  Effect  of  increasing  the  burner  wall  thermal 
conductivity 

In  order  to  achieve  uniform  temperature  profiles,  two  copper 
thermal  spreaders  were  installed  on  both  sides  of  the  burner  to 
increase  the  effective  axial  thermal  conductivity  [9] .  Each  copper 
spreader  was  5  cm  long  x  1  cm  wide  x  1 .8  cm  thick  and  weighed 
about  3.5  g.  Fig.  9  shows  the  steady  state  temperature  profiles 
with  the  copper  thermal  spreaders  installed.  The  higher  effective 
axial  thermal  conductivity  resulted  in  a  uniform  temperature 
profile. 

The  copper  thermal  spreaders  also  slowed  down  the  system 
start-up,  as  shown  in  Fig.  10.  The  slower  system  start-up  is  due  to 
the  increase  in  thermal  conductivity  and  mass  of  the  system.  The 
increased  mass  of  the  system  has  a  direct  effect  on  the  time  con¬ 
stant  for  start-up,  as  seen  from  Eq.  (1).  The  thermal  conductivity 
on  the  other  hand,  affects  the  time  constant  indirectly  by  mod¬ 
ifying  the  temperature  profile,  and  in  turn  the  local  rate  of  heat 
release.  As  methanol  and  air  enter  the  burner,  the  heat  released  by 
combustion  is  absorbed  by  the  burner  wall  and  causes  a  localized 
temperature  rise.  The  burner  having  low  axial  thermal  conduc¬ 
tivity  (without  copper  thermal  spreaders)  conducts  heat  slowly 
in  the  axial  direction,  and  therefore,  the  reaction  rate  at  the  hot 
spot  increases  faster,  leading  to  faster  system  start-up.  On  the 
other  hand,  the  burner  with  a  higher  effective  axial  thermal  con¬ 


ductivity  conducts  heat  faster  in  the  axial  direction  and  in  turn 
the  reaction  rate  is  slower  due  to  lower  temperatures  [19].  The 
slower  reaction  rate  and  the  increased  mass  result  in  a  slower 
system  start-up. 

3.5.  Microburner  integration  with  a  thermoelectric  device 

The  thermoelectric  module  used  was  Model  HZ-2  from 
Hi-Z  Technology,  Inc.  (San  Diego,  CA).  The  device  is  approx¬ 
imately  2.9  cm  wide  x  2.9  cm  long  x  0.5  cm  thick,  and  weighs 
13.5  g.  Two  pieces  of  1mm  thick  aluminum  nitride  electrical 
insulating  wafers  were  placed  between  the  secondary  thermal 
spreader  and  the  thermoelectric  device  to  prevent  electrical 
shorting  between  thermoelectric  elements.  To  ensure  good  ther¬ 
mal  contact,  thermal  grease  obtained  from  Hi-Z  Technologies 
was  placed  between  the  thermal  spreaders,  electrical  insulating 
wafers,  and  thermoelectric  module.  A  finned  copper  heat  sink, 
CoolWaveCPUCooler  from  Spire,  was  placed  on  top  of  the  ther¬ 
moelectric  module  to  enhance  convective  cooling  efficiency.  To 
further  ensure  thermal  contact,  the  entire  microburner/thermal 
spreader/thermoelectric/heat  sink  stack  was  placed  under  com- 
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Fig.  11.  Thermoelectric  transient  temperature  difference  and  power  vs.  time. 
Air  flow  rate  =  2000  seem,  9%  methanol  in  air,  AS  =  0.032  V  K_1 . 


pressive  force  using  a  conventional  C-clamp.  It  was  shown  in 
previous  work  that  this  compressive  force  has  a  positive  effect 
on  the  power  generated  from  the  system  [9]. 

Fig.  11  shows  the  temperature  difference  between  the  hot 
and  cold  sides  of  the  thermoelectric  as  a  function  of  time.  Also 
shown  is  the  calculated  power  using 

(  ASAT  \2 

P  =  *loadUioad  +  /W  (3> 

where  AS  is  the  difference  between  the  Seebeck  coefficients  of 
the  two  materials  of  the  thermoelectric  and  ATis  the  temperature 
difference  between  the  hot  and  cold  sides.  The  load  resistance 
Onload)  was  matched  to  the  thermoelectric’s  internal  resistance 
(/?int  =  5.3  Q)  to  calculate  the  maximum  power.  The  open  cir¬ 
cuit  voltage  of  the  thermoelectric  was  measured  at  steady  state 
under  different  air  flow  rates  and  methanol  mole  fractions.  For 
the  different  conditions,  AS  was  calculated  from  the  following 
equation  V=  AS  x  AT,  where  V is  the  open  circuit  voltage,  and 
we  used  the  average  AS  (0.032  VK-1)  to  calculate  the  power 
output  in  Fig.  11. 

As  shown  in  Fig.  11,  the  temperature  difference  between 
the  hot  and  cold  sides  increases  with  time  until  it  reaches  a 
maximum  around  6  min  where  the  output  power  is  maximum 
(0.65  W).  Then  the  temperature  difference  and  output  power 
start  decreasing,  and  they  reach  steady  state  around  15  min.  The 
reason  for  the  maximum  in  temperature  difference,  and  there¬ 
fore  in  power,  is  the  slow  heat  removal  from  the  cold  side.  As 
the  system  starts,  the  burner  starts  heating  up  the  hot  side  of  the 
thermoelectric  and  its  temperature  increases.  However,  heat  is 
also  conducted  to  the  cold  side  of  the  thermoelectric  where  it  is 
dissipated  and  removed  by  the  heat  sink.  This  creates  a  lag  time 
between  the  cold  side  transient  temperature  profile  relative  to 
the  hot  side  (data  not  shown).  Hence,  as  the  hot  side  temperature 
reaches  steady  state,  the  cold  side  temperature  keeps  increasing, 
resulting  in  a  maximum  in  the  temperature  difference  and  output 
power.  The  efficiency  was  calculated  using 

P 

efficiency  (%)  =  — x  100  (4) 

m  AH 


where  Pgen  is  the  maximum  power  generated,  m  is  the  mass  flow 
rate  of  methanol,  and  AH  is  the  heat  of  combustion  for  methanol 
(lower  heating  value).  The  maximum  efficiency  estimated  from 
the  data  was  1.1%. 

The  steady  state,  generated  power  for  the  methanol  system 
was  slightly  higher  than  that  of  previous  work  with  hydrogen 
and  propane  [9].  In  Fig.  1 1,  the  system  generates  approximately 
450  mW  with  a  temperature  differential  of  about  100  °C.  In  con¬ 
trast,  the  highest  power  output  from  the  previous  study  was 
approximately  420  mW  at  a  higher  temperature  differential  of 
1 10  °C.  This  is  most  likely  due  to  the  fact  that  the  thermoelectrics 
were  observed  to  have  different  Seebeck  coefficient  values,  as 
the  module  used  in  this  study  had  a  lower  operating  time.  Also, 
more  compressive  force  may  have  been  applied  to  the  methanol 
system,  which  was  shown  to  substantially  increase  the  heat  trans¬ 
fer  to  the  thermoelectric,  increasing  the  apparent  difference  in 
Seebeck  coefficients. 

With  respect  to  other  integrated,  portable  burner/thermo¬ 
electric  systems  utilizing  liquid  or  easy  to  liquefy  fuels,  Yoshida 
et  al.  has  demonstrated  power  generation  using  butane  with  a 
Laval  nozzle  ejection  system  [20].  However,  they  required  the 
use  of  resistive  heaters  to  start-up  the  system  and  were  unable  to 
sustain  catalytic  butane  combustion  autothermally,  i.e.,  without 
an  external  heat  supply. 

The  efficiency  of  the  integrated  microburner/thermoelectric 
system  still  needs  to  be  improved  in  order  to  become  a  compet¬ 
itive  alternative  to  batteries  in  low  power  portable  applications. 
However,  the  integrated  system  was  not  optimized.  Yet,  its  major 
advantage  is  its  simplicity  as  it  integrates  only  a  few  parts. 
Some  strategies  for  increasing  its  thermal  efficiency  include: 
(1)  use  of  a  second  thermoelectric  (e.g.,  on  the  other  side  of 
the  microburner)  to  utilize  more  effectively  the  heat  generated 
from  the  microbumer.  (2)  Heat  recovery  from  the  exhaust.  (3) 
Heat  removal  from  the  cold  side  of  the  thermoelectric  by  the 
reactants  through  a  micro-heat  exchanger.  This  approach  will 
pre-heat  the  reactants  (a  desirable  phenomenon)  while  increas¬ 
ing  the  temperature  difference  between  the  hold  and  cold  sides 
of  the  thermoelectric. 

4.  Conclusions 

Catalytic  self-ignition  of  methanol  over  Pt  nanoparticles  dis¬ 
persed  on  an  anodized  alumina  wafer  in  a  microburner  was 
demonstrated.  The  microburner  was  integrated  with  a  thermo¬ 
electric  device  and  the  system  was  also  shown  to  be  self- starting 
from  room  temperature  at  all  air  flow  rates  and  air/methanol 
ratios  studied.  The  system  start-up  time  was  highly  affected  by 
the  air  flow  rate  and  the  thermal  mass  of  the  system  but  nearly 
independent  of  the  mixture  composition.  The  system  reached 
steady  state  faster  at  higher  air  flow  rates.  These  findings  should 
be  fairly  generic  as  simple  theoretical  analysis  and  CFD  sim¬ 
ulations  indicate.  It  is  thus  proposed  that  high  flow  rates  and 
relatively  lower  fuel  content  should  be  used  in  order  to  reach 
steady  state  faster  with  minimal  fuel  utilization.  The  maxi¬ 
mum  power  generated  with  the  thermoelectric  was  0.65  W  and 
the  maximum  thermal  efficiency  was  1.1%.  This  is  the  high¬ 
est  reported  integrated  thermoelectric/microreactor  efficiency 
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reported  using  a  liquid  fuel.  The  efficiency  we  report  is  also 
close  to  that  of  an  integrated  fuel  reformer,  combustor,  fuel 
cell  system  [4].  Considering  the  simplicity  of  the  thermoelec¬ 
tric/microreactor  system,  it  is  probably  a  better  alternative  for 
low  power  (~1  W)  production  for  portable  devices. 
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Appendix  A.  Computational  fluid  dynamics  (CFD) 
simulations  for  the  effect  of  flow  rate  and  mixture 
composition  on  start-up  time 

In  order  to  verify  the  experimental  findings  and  assess  the 
insights  gained  from  the  simple  theoretical  analysis  presented 
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Fig.  A.l.  Schematic  of  computational  fluid  dynamics  (CFD)  geometry  and 
boundary  conditions  used. 


above,  CFD  simulations  were  carried  out.  The  microreactor  is 
modeled  as  two  parallel  plates  of  infinite  width  and  793  p,m 
(1/32")  thickness,  as  done  in  previous  work  [8].  The  combustion 
chamber  has  a  gap  size  of  400  jutm  and  is  1  cm  long.  Premixed, 
methanol/air  mixtures  are  sent  to  the  reactor  with  a  uniform  flow 
velocity  at  ambient  temperature  (300  K)  into  the  combustion 
chamber. 

The  CFD  package  Fluent  6.2  is  used  to  solve  the  spatially 
2D  time  dependent  continuity,  momentum,  species,  and  energy 
equations  using  a  finite  volume  method.  The  mesh  used  to 
perform  these  calculations  was  generated  in  the  commercial 
package  Gambit  and  contained  101  and  100  uniformly  spaced 


Fig.  A.2.  Maximum  wall  temperature  vs.  time  for  varying  composition  (a)  and  flow  rate  (c).  The  corresponding  normalized  transient  graphs  are  labeled  as  (b)  and 
(d),  respectively.  In  all  cases,  the  wall  thermal  conductivity  is  30  Wm-1  K-1  and  the  heat  loss  coefficient  is  20  Wm-2  K-1 
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nodes  in  the  axial  and  transverse  directions,  respectively.  The 
simulations  were  performed  on  2.4  GHz  Pentium  Xeon  proces¬ 
sors  of  2  GB  of  RAM  each. 

A  detailed  description  of  the  boundary  conditions  and  cal¬ 
culation  of  transport  properties  can  be  found  in  previous  work 
[21].  The  exterior  walls  were  modeled  assuming  Newton’s  law 
of  cooling,  as  seen  in  Fig.  A.l.  A  first-order  kinetics  model 
was  employed  along  the  catalytic  wall.  Given  the  lack  of  a  reli¬ 
able  kinetic  rate  expression  and  the  question  of  how  generic 
is  the  behavior  seen  experimentally,  a  pre-exponential  of  10 13 
1  s-1  and  zero  activation  energy  have  been  assumed.  The  lat¬ 
ter  choice  gives  self-ignition  in  the  transient  simulation.  No 
attempt  was  made  to  fit  the  experimentally  measured  values  (e.g., 
temperatures  and  times)  by  adjusting  the  kinetic  and  operating 
parameters. 

The  maximum  temperature  along  the  wall  versus  time  was 
monitored  for  varying  composition  and  flow  rate  to  validate 
the  assumptions  made  in  deriving  Eq.  (1).  Fig.  A.2a  and  b 
shows  the  effect  of  composition  at  constant  flow  rate  (gtot  =  1 
SLPM).  As  richer  fuel  content  mixtures  are  used,  the  heat 
released  from  the  catalytic  combustion  increases  and  thus 
both  the  steady  state  temperatures  and  heating  rates  increase. 
However,  when  the  temperatures  are  normalized  with  the  cor¬ 
responding  steady  state  temperatures,  the  transient  trajectories 
collapse  on  each  other.  Fig.  A.2c  and  d  shows  the  effect  of 
flow  rate  at  fixed  composition  (YMeOH  =  10%).  As  expected, 
the  steady  state  temperature  increases  with  increasing  flow 
rate.  In  addition,  the  system  reaches  steady  state  faster  at 
higher  flow  rates.  These  CFD  results  are  qualitatively  simi¬ 
lar  to  the  experimental  observations,  where  only  change  in 
flow  rate  had  a  significant  effect  on  the  time  to  reach  steady 
state. 
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